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We compare the single inclusive spectra of isolated photons measured at high transverse energy 

in proton-antiproton collisions at ^/s = 1.96 TeV with next-to- leading order perturbative QCD 

predictions with various parametrizations of the parton distribution functions (PDFs). Within the 

experimental and theoretical uncertainties, the Tevatron data can be reproduced equally well by the 

recent CTEQ6.6, MSTW08 and NNPDF1.2 PDF sets. We present also the predictions for isolated 

^1^' 7 spectra in proton-proton collisions at -^s = 14 TeV for central {y = 0) and forward {y = 4) 

D , rapidities relevant for LHC experiments. Different proton PDFs result in maximum variations of 

C/j . order ±30% in the expected _E^ -differential isolated 7 cross sections. The inclusion of the isolated 

r~^^ ' photon data in global PDF fits will place extra independent constraints on the gluon density. 
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I. INTRODUCTION 



The production of photons with large transverse energy {E^ ^ Aqcd ~ 0.2 GeV) in high-energy 

hadronic collisions is an interesting process in itself as a testing ground of the perturbative regime of 

Cs| ' Quantum Chromodynamics (pQCD) [l|, as well as a possible background for important new physics 

^ \ searches (such as e.g. often exemplified in the Higgs boson decay into two high-i?^ photons [2])- From a 

0^ , pQCD perspective, prompt photons issuing from hard parton-parton scatterings constitute a particularly 

clean testbed of perturbation theory in the collinear- [3|, |4| and kx- [5| factorization (or colour-dipole [g, Q ) 

approaches as well as of various logarithmic resummation techniques [8|. In addition, high-E'j, photons 

also yield valuable information about non-perturbative objects such as the parton distribution functions 

(PDFs) in the proton [3, ll^ and the parton-to-photon fragmentation functions (FFs) [ll|, II4I . At lowest 

order in the electromagnetic and strong coupling constants 0{aas), three partonic mechanisms produce 

prompt photons in hadronic collisions: (i) quark-gluon Compton scattering qg — >■ 7(7, (ii) quark-antiquark 

annihilation qq —^ 73, and (iii) the collinear fragmentation of a final-state parton into a photon^, e.g. 

qq — )■ qq — )■ 7X. The photons produced in the two first point-like processes are called "direct", the 

latter "fragmentation" photons. The Compton channel is particularly interesting as it provides direct 

^ ' information on the proton gluon distribution, g(x,Q^), which is otherwise only indirectly constrained via 

5^ I the derivative of the proton structure function dF2{x,Q'^)/dlogQ'^ ("scaling violations") in deep-inelastic 

scattering (DIS) e-p collisions [14|. 

Experimentally - see e.g. |15l . |l6[ for data compilations - many measurements of high-i?^ photon 
production have been carried out in the last 30 years in proton-proton (p-p) and proton-antiproton (p-p) 
collisions at fixed-target and collider energies spanning center-of-mass (cm.) energies of yG sa 20 - 
2000 GeV. Apart from results at fixed-target energies from the E706 experiment [T^, the agreement 
between data and theory is very good over nine orders of magnitude in the cross sections [16| . In order 
to identify the prompt photon signal out of the overwhelming background of photons from the decays 
of 7r° and t] mesons produced in the fragmentation of jets, one often requires the photon candidates to 
be isolated from any hadronic activity within a given distance around its direction. The corresponding 
measurements are then dubbed "isolated photons". The most recent^ iS^ -differential spectra for isolated 
7 - measured at RHIC {p-p at y/s = 200 GeV) ^ and Tevatron {p-p at V^ = 1-96 TeV) [H H^l - can 
be well reproduced by theoretical predictions based on Next-to-Leading Order (NLO) pQCD calculations. 



^ Though at first sight such a process is 0{a a^), the fragmentation function of the outgoing parton into a photon scales as 

a/ag and the corresponding cross section is of the same order as the direct process Il3|l . 
^ Not included in the aforementioned data compilations. 



Despite the existence of a few hundred data points, the inclusive and isolated photon data measured in 
p-p and p-p collisions have not been considered in the last ten years in the global-fit analyses performed 
to determine the proton PDFs. The last PDF parametrization to use the photon data was MRST99 [2l| . 
The theoretical difiiculties to reproduce the fixed-target E706 results at ^/s « 30 GeV [1^ motivated 
their removal from the PDF studies. Instead, Tevatron jet data were preferred over photon data in 
global analysis to constrain the (high-x) gluon PDF. The available comparisons of the latest measured 
isolated photon spectra to NLO have been carried out with just one or two PDF parametrizations, e.g. 
CTEQ6.1M 22] and MRST04 [23,], obtained a few years ago. Recently, new PDF sets have become 
available^- CTEQ6.6 [H], MSTW08 [H] and NNPDF1.2 ^ - which include the most up-to-date fits 
to data from deep-inelastic lepton-proton scattering and hadronic collisions as well as various theoretical 
improvements. The purpose of this paper is twofold. Firstly, we want to compare the existing Tevatron 
data to the theoretical predictions obtained using these latest PDF sets. Secondly, we want to study the 
expected sensitivity of near-future isolated photon spectra at the LHC to the same three sets of PDF 
parametrizations . 

The paper is organized as follows. In Section |ll] we succinctly remind the theoretical framework of 
our study and present the next-to-leading-order pQCD Monte Carlo (MC) programme jetphox that we 
have used. In Section IIIIl we present a comparison of the existing Tevatron isolated photon spectra with 
the JETPHOX results for various PDFs, theoretical scales and FFs. We then present in Section IIVI the 
predictions for LHC energies, at central rapidities (y = 0) accessible to the ALICE [2^, ATLAS [2^ and 
CMS ^Sfl] experiments as well as at the forward rapidities (up to y = 5) covered by LHCb [31[ . We discuss 
the associated uncertainties in the spectra and the sensitivity to the underlying PDFs. We summarize our 
main findings in Section FVl 

II. THEORETICAL FRAMEWORK 

A. Inclusive prompt photon production 

As aforementioned, two types of processes contribute to the prompt photon production cross section in p- 
p and p-p collisions: the 'direct' contribution, where the photon is emitted directly from a pointlike coupling 
to the hard parton-parton vertex, and the 'fragmentation' (sometimes also called 'anomalous') contribution, 
in which the photon originates from the collinear fragmentation of a final-state parton. Schematically, the 
differential photon cross section in transverse energy E^ and rapidity y can be written as [16J: 

da = da^.^ + da^^^^ ^ X! / ^^adxt Fa(xa; M^)Fb(a;6; /i^) x (1) 

a.,b=q,q,g 

do-2^{P'(,Xa,Xb;n^,fip,Hfj,)+ Y^ / —da^^^{p^,Xa,Xb,z;n^,fi^,fi^^.)D2{z;n'^^.^) 



c=q,q,g 



where Fa{xa, iJ-p) is the parton distribution function of parton species a inside the incoming hadrons h 
at momentum fraction Xa] Dj/i^(z, fj,,.) is the fragmentation function of parton fc to a photon carrying a 
fraction z of the parent parton energy (integrated from z^^n = XTCOshy, with xx = 2E^/y/s, to 1); and 
the arbitrary parameters /^^, /x^ and fi., are respectively the renormalisation, initial-state factorisation, 
and fragmentation scales which, loosely speaking, encode any residual dependence of the cross sections 
to higher-order contributions missing in the calculation. The study provided in this article relies on the 
calculation of both da^.^ and da^^^^ at next-to-leading order (NLO) accuracy [32] in the strong coupling 
as{^J,j^), i.e. all diagrams up to the order ©(aa^) are included, defined in the MS renormalisation scheme. 



3 During the completion of this work, the NNPDF collaboration released new PDF sets (NNPDF2.0) [l^. Although the 
gluon PDF has changed, it is consistent with the NNPDF1.2 version used here, and we do not expect significant difli'erences 
on the photon production predictions. 



The results of the NLO calculation of da^.^ have been known for a long time 3] . The calculation of the NLO 
corrections to da^^^^ became also available later [J, [ij, [S^] ■ We note that the distinction between da^-^ and 
dcT.^ is arbitrary (only its sum is physically observable): a typical case is a bremsstrahlung process from 
a final-state quark which, depending on the scale, can be considered as "fragmentation" or as "NLO direct" . 

More recently, expressions involving the resummation - at next-to-leading (NLL) or even next-to- 
next-to-leading (NNLL) logarithmic accuracy ~ of threshold and recoil contributions due to soft gluon 
emission, which are large close to the phase space boundary where the E^ of the photon is about half of 
the center-of-mass energy {xt ~ "^E^j %fs — >■ 1) have been obtained for da^-^ and da ^^^^ iS|. The effect 
of this resummation is important at very large^ photon Erj, extending down to values of xt ~ 10~^, and 
provide a much reduced /i^, and /i^ scale dependence than the NLO approximation. In any case, since we 
are interested in the low-x region of the PDFs, we do not consider such effects, which are not implemented 
in JETPHOX, in this study. 

In Figure [1] we show the kinematical region covered by the photon measurements in the (xt, E'^') plane 
equivalent to the DIS (x^Q"^) plane for hard scattering hadronic processes. The quoted x values are typical 
XT values for measurements at midrapidity and x k. xt • exp(-2/) at forward rapidities. At midrapidities, 
photon production at the LHC will probe values 20 times smaller than at Tevatron because of the factor ^7 
increase in the cm. energy and the lower E^ values reachable (e.g. down to E^ « 5 GeV at ALICE). At the 
forward rapidities and low E^ k, 5 GeV covered e.g. by LHCb, one can probe the gluon distribution down 
to very small x « 10^^. The combined photon data at all energies - about 350 data-points collected so far 
plus O (100) new data-points expected from the four experiments at the LHC - provide thus constraints 
of the PDFs in a wide (x, Q^) range. 
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FIG. 1: Kinematical region probed by existing prompt photon measurements at fixed-target (Fermilab) and collider 
(ISR, RHIC, SppS, Tevatron) energies, and expected range probed at the LHC at central (j/ = 0) and forward 
(j/ = 2 - 5) rapidities. 

In Figure [5] we show the relative contributions of each one of the three subprocesses (Compton, 



* Note that at the highest momenta probed, there are also additional corrections due to electroweak boson exchanges, which 
decrease the photon yields by about 10-20% within E' ■^ \ - 2 TeV at the LHC [H. 



annihilation and fragmentation) to prompt photon production at Tevatron (y = 0) and LHC (y = 0, 
and y = 4) as a function of the photon i?^ . [Again, the plots are to be taken indicatively as the weights 
of the different components are scale-dependent.] They have been obtained selecting the correspond- 
ing Feynman diagrams at NLO with the jetphox MC, setting all scales to jj, — E^, and using the 
CTEQ6.6 parton densities and the BFG-II parton-to-photon FFs for the fragmentation photons. In 
the low Ej, region of the mid-rapidity spectra (below 15 GeV at Tevatron, and 40 GeV at the LHC) 
the fragmentation component dominates the cross sections. At the Fermilab collider, quark-gluon 
scattering is the dominant component up to E^ « 120 GeV, beyond which the annihilation of valence 
(anti)quarks from the (anti)proton beams play a preeminent role. At the LHC, the Compton process 
dominates for all E.^'s above 45 GeV at y = 0. At forward rapidities at the LHC, between 50% to 80% 
of inclusive prompt photons are produced from the coUinear fragmentation of a parton for all relevant E^ 's. 
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FIG. 2: Relative contributions of the quark-gluon Compton, qq annihilation and fragmentation subprocesses in 
NLO prompt photon production at the Tevatron (left), LHC midrapidity (center) and LHC forward rapidity (right) 
obtained with jetphox (CTEQ6.6 PDF, ^ = E^^, and BFG-II FFs). 



B. Isolated photon production 



In order to strongly suppress the overwhelming background of secondary photons coming from the 
electromagnetic hadron decays (mainly tt'' and 77 mesons) produced inside partonic jets, collider experiments 
require often an isolation criterion on the photon candidates. A standard isolation requirement is that 
within a cone about the direction of the photon defined in rapidity y and azimuthal angle by 



R 



\l {y - Vif + {<i} ~ <l^if 



(2) 



the accompanying hadronic transverse energy E]^ is less than a given value 



rphad ^ rpmax 



(3) 



R is usually taken between 0.4 - 0.7 and E"^""^ is specified either as a fixed value, or as a fixed fraction 
£h (e.g. often 10%) of the photon E^ . Cross sections for producing isolated photons have been proven 
to fulfill the factorization theorem, and are finite to all orders in perturbation theory for non zero R and 

^max [3l|_ 



If one ignores the "underlying event" activity in the collision - due to soft interactions of spectator 
partons and beam remnants [35| - the leading-order qf^-Compton and gg-annihilation photons are emitted 
in a region completely free from hadronic activity (the recoiling quark or gluon is emitted at vr from it). 



On the other hand, the fragmentation component is usuahy accompanied by other particles issuing from 
the hadronization of the same parent parton (jet). Yet, the isolated cross section measured experimentally 
cannot be automatically identified with the direct cross section calculated at the Born level. First, higher 
order terms originating in the non-collinear fragmentation of partons also contribute to the isolated cross 
sections. Second and most important, although isolation cuts such as Eqs. ([2]), ([3]) reduce the da^^^^ 
contribution, a fraction of fragmentation photons with z > 1/(1 + Sh) survive the cuts [16[. The average 
2;-value for fragmentation photons is (z) < 0.6 at the LHC and (z) « 0.7 at the Tevatron [ij], and a 
typical isolation energy cut of Sh — 0.1, corresponding to 1/(1 + Sh) > 0.9, suppresses about 60 - 80% of 
d'^fr-ag (thls valuc Is E^- and scale-dependent). We see this in more detail in Fig. [3] where we show the 
subprocesses contributions to the isolated photon cross section. At variance with Fig. [2] for the inclusive 
prompt-7 case, we can see that a very significant part of the fragmentation component is suppressed after 
applying typical isolation cuts {R = 0.4, Eh = 0.1). The Compton process now clearly dominates the 
photon yield below E^ « 120 GeV at Tevatron and accounts for about 3/4 of isolated 7 production for 
all £'j,'s at the LHC. These results confirm the interesting potential of isolated photon hadroproduction 
to constrain the gluon density. 
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FIG. 3: Relative contributions of the quark-gluon Compton, qq annihilation and fragmentation subprocesses in 
NLO isoZaied photon production at the Tevatron (left), LHC midrapidity (center) and LHC forward rapidity (right) 
obtained with jetphox (CTEQ6.6 PDF, fi = E'^, BFG-II FFs) for an isolation radius R = 0.4 and a hadron fraction 
of the photon energy of e^ = 0.1 inside the cone. 



JETPHOX Monte Carlo 



The present study relies on the implementation of the photon NLO calculation of both dcr^^^ and 



dff^^^g in the jetphox Monte Carlo (MC) programme [32 



The main advantage of the JETPHOX MC 



is that one can easily account for any kind of experimental cuts (e.g. on kinematics and/or isolation) 
implementable at the partonic level. In addition, one can match naturally the binning of experimental 
data by histogramming of the partonic configurations generated. All the NLO results provided in the 
following sections are obtained using^ as{Mz) = 0.118, with up to 5 active quark flavours (the lowest E^ 
considered in this work is close to nii, « 4 GeV/c^). We switched off the box diagram .gg — >■ 37 in the 
calculations because its contribution to the single inclusive spectrum is found to be of just a few percent. 
The CTEQ6.6, MSTW08 and NNPDF1.2 PDFs were interfaced to jetphox via the lhapdf (version 5.7.1) 



5 Note that although each PDF set uses a slightly diflferent reference value: as{Mz) = 0.118 (CTEQ6.6), as{Mz) = 0.119 
(NNPDF1.2), and as(Mz) = 0.12018 (MSTW08), the resulting differences in the 7 cross sections obtained using the slightly 
different coupling choices are very small. 



package [36|. Whenever the scales /i^, ^p and /i , , , are given a common value, the latter is noted /i hereafter. 

In JETPHOX the NLO calculations are performed at the parton level and do not account for hadronisation 
effects which can be potentially large at low E^ . In addition, they do not include the soft hadronic activity 
from the underlying event (parton spectator collisions, beam remnants ...) whose transverse energy can 
also fall inside the isolation cone. However, the experimental data at Tevatron have been corrected for such 
effects. In the case of the CDF analysis, the correction is found to amount about 9% of the photon yield, 
independent of its E^ above 30 GeV [T^l ■ Studies at the LHC using various tunes of the underlying event 
with PYTHIA [33] indicate that the loss of isolated photons due to hadronic activity within the isolation 
cone is of about 20% at E^ « 5 GeV and negligible above 70 GeV (for radius R — OA and hadronic energy 
fraction £,, = 0.1) ^. 

III. COMPARISON OF TEVATRON DATA TO NLO PQCD 

During the Tevatron Run-II, the D0 and CDF collaborations collected respectively 0.32 fb~^ and 2.5 fb~^ 
worth of photon data in p-p collisions at y^ =1.96 GeV [la, (20]. Such data sample allowed both experi- 
ments to measure isolated prompt photons in the E^ range from roughly 20 GeV to about 400 GeV. The 
detailed conditions of the measurements are as follows 

• D0: Kinematics: E'J = 20 - 300 GeV, \y''\ < 0.9. Isolation criterion: R = 0.4, en = El<"^ / E^ < 0.1. 

• CDF: Kinematics: £'J = 30 - 400 GeV, jy^l < 1.0. Isolation criterion: R = 0.4, E'^""^ < 2 GeV. 

We have used the kinematics and isolation cuts in jetphox corresponding to each measurement and 
calculated the resulting isolated photon spectrum with the same E^ binning as the experimental spectra. 
The calculations have been repeated for three PDF sets (CTEQ6.6M3ISTWO8, NNPDF1.2), three scales 
{^i=0.5EI, eI, 2EI) and two photon FFs (BFG-I and BFG-II) jll|. The results are presented and 
discussed in the next subsections. 



A. Parton distribution functions dependence 

The D0 and CDF isolated photon spectra are compared in Fig. 2] to the corresponding jetphox pre- 
dictions for three different PDF parametrizations (CTEQ6.6, MSTW08 and NNPDF1.2). The theoretical 
scales have been set to ji — E^ and the parton-to-photon FFs to the BFG-II set. The three theoretical 
spectra are very similar and they all agree well with the data within the experimental uncertainties given 
by the error bars (resp. bands) representing the statistical (resp. systematical) errors. The isolated photon 
spectrum obtained using MSTW08, for both sets of experimental cuts, is mostly in between the spectra 
obtained with the two other PDFs, except at small- £■ where MSTW08 results in a somewhat lower cross 
section, closer to the data. This is better seen in Fig.Elwhere we plot the ratio data/theory calculated with 
the three choices of PDF. On average, the NNPDF1.2 (resp. CTEQ6.6) prediction is a few percents higher 
(resp. lower) than the MSTW08 distribution. As aforementioned, all PDFs reproduce well in magnitude 
and shape the experimental spectrum within the experimental uncertainties, dominated by the energy-scale 
uncertainty (orange bands in the plots). 

A simple x^ analysis of the full spectrum to the theoretical predictions, accounting for the experimental 
uncertainties (statistical and systematical added in quadrature) yields y^ /ndf below 1 for the three PDF 
predictions. For the default scale-choice shown, ^ — E , MSTW08 reproduces a bit better both the CDF 
and D0 results. One has to note, however, that in the E range below ~40 GeV a shape difference is 
apparent, with all NLO curves systematically under-predicting the steepness of the cross section by up to 
around 40% in the lowest E^ bin. The source of such a disagreement, which is about a 2a effect in the case 
of D0, is unclear. Theoretically, the inclusion of small-x resummation effects [33] does not help to improve 
the data-NLO agreement. 
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FIG. 4: Isolated photon iJ^ -differential cross section at mid-rapidity in p-p at -^s — 1.96 TeV measured by D0 
(left) and CDF (right) compared to jetphox {fi = eI, BFG-II) for three different PDFs: CTEQ6.6, MSTW08 and 
NNPDF1.2. The (orange) bands indicate the experimental systematic errors. 
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FIG. 5: Ratio of isolated photon midrapidity spectra in p-p at ^/s = 1.96 TeV measured by D0 (left) and CDF (right) 
over the jetphox {fj. = eI, BFG-II) predictions for three PDFs: CTEQ6.6, MSTW08 and NNPDF1.2. For clarity, 
only the MSTW08 ratio is shown with propagated statistical errors from the data. The lines between data/NLO 
points are to guide the eye. The (orange) bands indicate the systematic uncertainties of the measurements. 



B. Scale dependence 



As mentioned in Section III Al the truncation of the as expansion of the pQCD cross section for prompt 
photon production, Eq. ([2]), at a given order introduces an arbitrary dependence on three unphysical 
scales: the renormalization scale /^„ which appears in the evolution of the coupling, the factorization scale 
fj,p associated to coUinear singularities in the PDFs, and the fragmentation scale fx^^ related to collinear 



8 

singularities in the FFs. Since those scales are unphysical, the theory can be considered reliable only where 
the predictions are stable with respect to variations of them. Usually scale variations of a factor of two 
around the "natural" physical scale of the scattering process, often taken as the E^ of the photon, are 
considered: /i — E^ /2 - "^.E^ . In Fig. [51 we show the ratio of the JETPHOX isolated 7 spectra at Tevatron 
obtained with scales set to ^ = E /2 and jjl — 2E over the same spectrum for /i = £■ . In all three 
cases, we use the CTEQ6.6 PDFs^ and the BFG-II FFs. The sensitivity to the changes in the theoretical 
scale /x is of about ±10% in the whole ^^.-range. The choice of a scale value fi = E^ agrees better with 
the experimental spectrum for the MSTW08 and CTEQ6.6 PDFs whereas the NNPDF1.2, whose central 
spectrum is a bit higher than the two others, favours a larger /i — 2E scale. 
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FIG. 6: Scale uncertainty of the jetphox NLO predictions for isolated photon production at Tevatron (CDF 
isolation cuts) plotted as the ratio of the spectrum obtained with jj, = 0.5E^ , 2E^ over that with /i — E^. 



C. Dependence on the parton-to-photon fragmentation functions (FFs) 



The production of fragmentation photons in hadronic collisions is completely encoded in the parton-to- 



photon FFs, D^/fc(2;,/ij:j,) obtained from fits [ll|,|40| to photon production data at LEP (e"* 



99(5) -> 



7X) [41[. In JETPHOX two choices of FFs are available: the BFG set I ("smaU gluon") and II ("large 
gluon") [il|. Although both FFs reproduce the LEP data, mostly sensitive to the high z = p-y/pparton 
range of the FFs, at small z the gluon-to-photon fragmentation function of set I is significantly lower than 
the one of set II. Otherwise, the quark-to-photon fragmentations are identical in BFG-I and BFG-II. We 
have run prompt and isolated photon production in the Tevatron kinematic range with jetphox at NLO 
with the two FF sets. The inclusive prompt 7 spectrum obtained with BFG-I is just a few percent smaller 
than the one obtained with BFG-II, below E^ « 80 GeV. The corresponding BFG-I and BFG-II isolated 
photon spectra are virtually identical. This is not unexpected since a significant fraction, up to 80% (see 
Fig. [3]), of the fragmentation photons are removed by standard isolation cuts R = 0.4, Sh = 0.1- 



The scale-dependence obtained with the two other PDF parametrizations are similar. 



IV. ISOLATED PHOTON PRODUCTION AT THE LHC 

At the top LHC center-of-mass energy of 14 TeV, the production cross section of isolated photons, <j{E^ > 
10 GcV)|.y=o = 300 nb at NLO, is seven times higher than at Tevatron, a{El > 10 GeV)|j^=o = 40 nb. 
This will allow one to carry out high-statistics measurements even with moderate integrated luminosities, 
without (or with minimal) pileup. The ALICE, ATLAS, CMS and LHCb experiments at the LHC have 
photon reconstruction capabilities with electromagnetic calorimetry in various rapidity ranges. The two 
general-purpose high-luminosity detectors ATLAS [2g| and CMS [30| cover'^ \y\ < 3 from E^ w 10 GeV up 
to very large E^ sa 1 TeV, whereas ALICE [2^ covers \y\ < 0.7 in a range of moderate E'^ « 5 - 100 GeV. 
LHCb [3l|, on the other hand, can measure photons within 2 < 77 < 5 albeit up to relatively moderate 
Er^ PS 20 GeV, as the calorimeter dynamic range is mostly optimised for low energy photons coming 
from i?-mesons radiative decays. As possible benchmark measurements, we will consider isolated (radius 
R = 0.4) photon production at y = (for ALICE, ATLAS and CMS) and at y == 4 (for LHCb*^). The 
JETPHOX predictions for central and forward rapidities for three different PDF parametrizations (CTEQ6.6, 
MSTW08 and NNPDF1.2) are shown in Fig. H 



> 10^ 


^ 




^ 10"^ 
Q. 

>. 10^ 
'h 10 


j pp^yi3<„X,\/s=14TeV,y=0 
i — CTEQ6.6 
J — MSTW08 
= \ — NNPDF1.2 

^\ JETPHOX 1.1 : n=EJ, FF=BFG-II 


1 


: \ Y isolation: R=0.4,ef,=0.1 


10-^ 


r \. 


10-2 


r \.^^ 


10-^ 


r ^^"^^^-.^^^^^^^ 


10-* 


r ^^'"^""^^ 


in-5 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



100 200 300 400 500 600 700 800 900 1000 

e; (GeV) 



> '^ 


: 


0) 



S10S 


PP^YisoiX,\/s=14TeV, y=4 

^ — CTEQ6.6 


>> 


^ — MSTW08 


LLI 

.^10' 


-\ — NNPDF1.2 


■D 


: \ 


10* 


r \v 


10^ 


r ^-\^^ 


10^ 


r JETPHOX 1.1:h=E|,FF=BFG-II 




Y isolation: R=0.4,eh=0.1 


10 





10 15 20 25 30 35 40 45 

E| (GeV) 



FIG. 7: Predictions of the isolated photon spectrum in p-p collisions at -^s = 14 TeV obtained with jetphox 
(scales fixed ai fi = eI and FFs set to BFG-II) for three different PDF parametrisations (CTEQ6.6, MSTW08 and 
NNPDF1.2) for central {y — 0, left) and forward (y = 4, right) rapidities. 

Differences of a few tens of percent in the photon spectra obtained with different PDFs are not 
noticeable in logarithmic plots such as those in Fig. [71 In Figure |5] we show the fractional differences, 
(1±PDF1/PDF2), between the jetphox predictions for isolated photon spectra {R = 0.4 and Eh = 0.1) 
obtained at y = in p-p collisions at ^/s = 14 TeV for the various PDF sets. They have been obtained 
with the CTEQ66, MSTW08 and NNPDF1.2 for fixed scales at /i = eI and BFG-II FF. In general, 
the agreement among different PDFs is very good except at very low and very high E^ . The choice 



''' We do not consider here the possibilities covered by the ATLAS and CMS forward hadron calorimeters up to \y\ < 5. 
** The chosen rapidity y = 4, lower than the maximum LHCb y = 5 value, accounts for the effect of the R = 0.4 isolation 
radius and possible additional fiducial cuts. 
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of PDF results in up to ^15% variations of the photon yields below E^ « 15 GcV (left plot). Similar 
PDF uncertainties are obtained within one single PDF parametrization. As an example in the right plot 
we show the envelope of maximum relative differences at each E between each one of the ten^ first 
NNPDF1.2 rephcas and the defauh (central) NNPDF1.2 set. 
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FIG. 8: Fractional differences between the isolated photon spectrum in p-p collisions at -^s = 14 TeV at central 
rapidities (y = 0) obtained with jetphox (scales fixed at fi = E^, FF set to BFG-II) using: (a) three different 
PDF parametrizations (CTEQ6.6, MSTW08 and NNPDF1.2), and (b) ten NNPDF1.2 replicas over the reference 
parametrization of the same PDF set. The dashed lines indicate a ±15% uncertainty. 

The expected PDF sensitivity for measurements of forward isolated photons {y = 4) in the transverse 
energy range E^ =5 — 50 GeV is shown in Figure [S] Variations in the photon yields of up to ~30% 
depending on the PDF set can be seen at the lowest E^ w 5 GeV (left plot). Intra-PDF uncertainties seem 
to be somewhat lower (around ~15% within 10 replicas of the NNPDF1.2 parametrization, right plot). 

As in the Tevatron case, the dominant source of experimental uncertainty (orange bands in Fig. 2]) 
in any (isolated or not) photon measurement at the LHC will likely come from the uncertainty in the 
absolute energy scale of the calorimeters. A typical 1.5% uncertainty in the photon energy calibration 
will propagate into about 10% uncertainties in the photon cross sections for a steeply-falling photon 
spectrum with inverse power-law exponent n ss 5 - 8. Altogether, ten to fifteen percent uncertainties in 
the final cross sections (not accounting for an overall luminosity normalization error) are to be expected 
at the LHC. Those values are usually better than the uncertainties linked to jet measurements where the 
energy-scale uncertainty, involving a good knowledge of the hadronic calorimeters calibration, is often 
poorer than for photons. 

In the theoretical side, as seen in the case of the Tevatron predictions fSection lllll) . the choice of the 
parton-to-photon FF plays no role on the final isolated photon cross sections and we are left with the 
choice of theoretical scales /z^, fi^ and /i,, as the only important source of uncertainty. In Figure [TUl 
we present the fractional differences, (1 ± /i'//i), between the predicted isolated photon spectra obtained 
with JETPHOX (MSTW08 PDFs and BFG-II FFs) with "default" theoretical scale {fi = ED and those 
obtained with /i' = E^ /2 - 2E^ . At low i?^ the uncertainty in the cross sections linked to scales variations 
amounts to up to ±20% but it saturates at around ±10% above E sa 15 GeV at central and forward 



^ Although this is a relatively small number of replicas, running .JETPHOX at NLO is a quite time-expensive operation from 
the computational point of view. 
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FIG. 9: Fractional differences between the isolated photon spectrum in p-p collisions at ^/s = 14 TeV at forward 
rapidities (y = 4) obtained with jetphox (scales fixed at p = E^, FF set to BFG-II) using: (a) three different 
PDF parametrizations (CTEQ6.6, MSTW08 and NNPDF1.2), and (b) ten NNPDF1.2 replicas over the reference 
parametrization of the same PDF set. The dashed lines indicate a ±10% uncertainty. 
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FIG. 10: Scale uncertainty plotted as the fractional differences between the isolated photon i?^ -differential cross 
section in p-p at y^ = 14 TeV at central {y — 0, left plot) and forward {y = 4, right plot) rapidities with scales 
fj, = E^ compared to that with scales: jj, = 0.5E^, 2E^. The dashed lines indicate a ±10% uncertainty. 



As aforementioned, global PDF fits do not reliably determine the gluon for x < 10^^ at low, yet 
perturbative, Q^ < 20 GeV^ scales [43] ■ This is due partly to the lack of precise structure function data 
for X < 10~^ and mainly due to the fact that the copious DIS F2 results probe the quark distribution, 
while the gluon density is constrained by the evolution, log(Q^) dependence, of these data. In the low x 
region the available Q^ range decreases since x and Q^ are correlated {x ex Q^ / s) and the accuracy of 
the gluon determination becomes worse. The main motivation to use the LHC isolated photon spectra 
to help constrain the low- a; gluon density in the proton is summarized in the left plots of Figs. [5] and [5] 
Different PDFs predict isolated photon yields which, at the lowest E values, differ by up to ±30%. Most 
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of such differences arise from the different small-x gluon densities implemented in the three PDF sets. 
Indeed, as seen in Fig. [S] the Compton quark-gluon scattering clearly dominates isolated 7 production at 
the LHC and, as indicated by Fig. [TJ the LHC data at forward (central) rapidities probes parton fractional 
momenta around x — 10^'* (resp. x — 10^'^). It is not unexpected that the largest PDF-dependent 
variations of the yields are observed in the low E^ range since it is indeed at moderate virtualities 
that the gluon density is more uncertain. The minimum photon E^ w 5 GeV considered here is pretty 
optimistic from an experimental point-of-view - the reconstruction and isolation of photons becomes more 
challenging with decreasing transverse energies - but it still corresponds to relatively large scales of order 
Q^ « 25 GeV^. Having the possibility to measure isolated photons with lower, but still perturbative, mo- 
menta around E^ « 3 GeV would provide even more stringent constraints on the low- a; gluon distributions. 

One can take full advantage of the measured isolated photon results to better determine g{x,Q^), by 
including them in future global-fit analyses. Indeed, as seen in Fig. I10[ the scale uncertainties at low 
E^ are relatively large (±20%) and thus a simple direct spectral comparison of the experimental E^- 
differential cross sections to the theoretical predictions will not unambiguously indicate the preferred PDF 
parametrization. The latest works of the NNPDF collaboration [23| use, for example, about 3500 exper- 
imental points to carry out their PDF fits, including DIS (F2, Fl and i^'charm^ ^^^^ proton- (anti)proton 
(Drell-Yan, vector-boson and jets) results. Yet, in the small- a; region the available experimental statistics 
is quite limited: just about 19 (126) data-points below x = 10""' (resp. x = 10"'^), among which only Fl is 
directly sensitive to the gluon. In these regions where there are little or no experimental constraints on the 
gluon, the most important source of the PDF uncertainty is due to the parametrization bias. Inclusion of 
new data-sets is basic to improve the central value of the PDFs and reduce their associated uncertainties. 
In that respect, including 10 ~ 20 data-points from forward isolated photon spectra at ^/s = 14 TeV, will 
place new direct experimental constraints on the gluon distribution at a; = 10""*. The effect of these data 
on future PDF parametrizations might be larger if the existing measurements have inconsistencies which 
previously have not been accounted for. All in all, this work indicates that the combined photon data at 
all energies - about 350 data-points collected so far plus O (100) new data-points expected from the four 
experiments at the LHC - provide new useful constraints of the gluon distribution in a wide {x, Q^) range. 

V. CONCLUSIONS 

We have compared the most recent measurements of isolated photon spectra in p-p collisions at 
\/s = 1.96 TeV [la, 12Q| with next-to-leading-order calculations with updated parton distributions 
functions (PDFs) in the proton. Previous data-theory comparisons [la, [T^ [20| used relatively older 
PDF sets. The i?^ -differential spectra at Tevatron can be equally well reproduced by the CTEQ6.6, 
MSTWQ8 and NNPDF1.2 global fit parametrizations within the existing experimental and theoretical 
uncertainties. Differences in the isolated photon cross sections for different PDF sets are in the range 
5% - 10%, with the MSTWQ8 spectrum being somewhat in between the spectra computed with the two 
other PDFs. Theoretical uncertainties linked to the scales choices /i^ = E^/2 — 2E^ result in around ±10% 
S^ -independent changes in the NLO cross sections. The choice of parton-to- photon fragmentation function 
(FF) does not introduce any additional theoretical uncertainty as the isolation criteria help to keep the 
contribution of fragmentation photons to the total prompt photon yield well below the fifteen percent 
at all transverse energies, and differences between FFs are a few percent of this remaining 15% contribution. 

We have also presented the NLO isolated photon spectra expected in p-p collisions at the top LHC 
energy for central {y — 0) and forward (y = 4) rapidities and determined the associated theoretical 
uncertainties. At \fs — 14 TeV, the prompt photon cross sections (above 10 GeV) are more than a factor 
of seven larger than at Tevatron, and the parton fractional momenta probed are relatively small, down 
to a; « 10~^ (resp. 10~^) at mid (resp. forward) rapidity. Since the quark-gluon Compton process is 
found to represent about two thirds of the total isolated photon cross section at the LHC for standard 
isolation cuts, such a measurement promises to provide an interesting direct measurement to the relatively 
unconstrained low- a: gluon density g{x, Q^) in the proton. The three PDFs studied yield central values 
of the photon spectrum which differ at most by ±15% at y = and by ±30% a,t y — 4 in the low-£^^ 
region of the spectra, whereas the experimental (associated with the calibration of the energy scale in the 
electromagnetic calorimeters) and the theoretical (linked to the choice of the renormalisation, factorisation. 
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and fragmentation scales) uncertainties are of O (10 — 20%). To finalize, we have presented arguments to 
motivate the inclusion of the LHC photon data in future PDF global-fit analyses. 

In summary, given (i) the high-quality isolated photon measurements available or expected at collider 
energies, (ii) the good agreement between the existing RHIC and Tevatron (except for the lowest E^ 
bins) data and NLO calculations, and (iii) the large statistics photon data expected to be collected at 
various rapidities at the LHC, it is worth to reconsider the inclusion of the combined photon collider data 
measured in proton- (anti)proton collisions - about 350 data-points collected so far plus O (100) new data- 
points expected from the four experiments at the LHC - in global-fit PDF analyses. Such data provide 
direct extra constraints on the gluon PDF in a wide (x, Q^) range. 
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